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Abstract 

Experimental data obtained by Lughi and Clarke [J. Am. Ceram. Soc., to be published] are 
compared to the theory describing disorder effects on optical phonons. Sharpening and vanishing 
of the asymmetry of the Raman lines after annealing are attributed to a decrease in short-range 
disorder. The parameters of disorder (such as the phonon-defect coupling and the correlation 
length) are defined in the comparison. 
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I. INTRODUCTION 



Raman scattering has been used for years to measure the center-zone phonon frequencies. 
But it was understood quite recently that studies of the Raman line shape can also provide 
widespread information relative to disorder, for instance, the isotopic compositions jjl], the 
stacking faults , the strain fluctuations near so on. The phonon 

scattering processes by defects of various types result in increasing the phonon relaxation 
rate and consequently the Raman line width. 

A novel field of Raman spectroscopy, the transformation kinetics, was put forward by 
Lughi and Clarke [5|. They investigated the phase transformations of yttria-stabilized zir- 
conia.Zirconia Zn0 2 has the wide range of applications, including traditional structural 
refractories, fuel cells, and electronic devices such as oxygen sensors. Because of these tech- 
nological implications, the characterization of zirconia is of particular interest. Zirconia has 
three zero-pressure polymorphs: the high-temperature cubic (c) phase (stable between 2570 
K and the melting temperature 2980 K), the tetragonal (t) phase (stable between 1400 and 
2570 K), and monoclinic (m) phase below 1400 K. 

The unit cell of the c phase contains one formula unit. There is only one triply degenerated 
Raman active mode at 640 cm -1 From the symmetry analysis, it is suggested 

that the c — t transformation of Zn02 results from a condensation of the Brillouin zone 
boundary phonons. Indeed, it was found in the model calculation of phonon curves, 

that the phonon frequency u;(X2~) of the c— phase is close to zero. The volume of the 
elementary ce 
active modes 

In the work |5|], zirconia coatings were obtained by electron beam evaporation and then 
annealed. The X-ray diffraction and Raman spectroscopy measurements were performed 
before and after annealing. They give a possibility to control the phase transformation in 
the samples. The obtained results are very interesting because they demonstrate the gigantic 
effect of annealing on both the width and asymmetry of the Raman line. 

In the present paper we compare the results jflj with the theory describing the disorder 
effect on the optical phonons. In Sec. II, a sketch of the theory is presented. The comparison 
of the theory with the experiment ^ are given in Sec. III. In Sec. IV, concluding remarks 
are summarized. 



1 is doubled in the c — t transformation. In the t phase, there are six Raman 
(their frequencies are given in Table). 
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II. THEORETICAL BACKGROUND: PHONON ELASTIC SCATTERING BY 
STATIC DISORDER 

The scattering of phonons by disorder (defects or strain fluctuations) results in the phonon 
width and shift (see hj for review ). To estimate the effect of disorder we can use the 

n 

quantum-mechanical formula [16J: 



giving the phonon width as a function of the frequency and wave vector. Here q is the 
concentration of defects, v(k— q) is the phonon-defect interaction, u is the phonon frequency 
in the initial state, Uq ± s 2 k 2 is the squared phonon frequency in the final state (the "+" 
sign if the phonon branch has a minimum at ujq and the "— " sign for the maximum), and 
8(x) is the Dirac ^-function expressing the conservation low in the phonon elastic scattering. 
The summation is performed over the 3-dimensional wave vector k in the case of phonon 
scattering by the point defects under consideration. 

Let us consider the short-range defects with the interaction having the form of the 8- 
function in the real space. Then, the potential v is independent of the wave vector, i.e. 
v (k — q) = vq with a constant value Vq, and the integration can be done explicitly. For the 
phonon minimum, the contribution of the scattering by defects appears only in the region 

UJ 2 > UJn'- 



In the case of the phonon branch maximum, we have the disorder contribution only if 

UJ 2 < Uq. 

Expression (J2J) has a transparent physical meaning: the phonon elastic scattering from 
defects produces some phonon width in that frequency range, where the appropriate final 
phonon states occur. For instance, in the case of the minimum of the phonon branch at 
uj = ujq, the final phonon states are only for uj > ujq- Of course, this effect manifests itself as 
a line asymmetry in the Raman light scattering on phonons. As a result, the high-frequency 
wing of the line drops more slowly (than the low-frequency one) for a phonon branch having 
the minimum at the center of the Brillouin zone. 

I should emphasize that the frequency dependence of the phonon-defect width is de- 
termined by the dimensionality of defects. As we see from Eq. (J2J), this is a square-root 




(1) 



A' 




(2) 
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dependence for point defects. In the case of linear defects such as dislocations (or the plane 
defects - crystallite boundaries or stacking faults), we have the 2d (or Id) integration in Eq. 
(HJ) which gives 



r(w) = T nat + -^\v \ 2 6(uj 2 - co 2 ), linedefect, (3) 



T{u) = T nat + -^\v \ 2 (u 2 - to 2 )- 1/2 , plane defect, (4) 

where 6(x) is the Heaviside step-function. 

The singularities in Eqs. (jSJ) - as well as the weak singularity in Eq. (j2J) arise, because 
the width of final states was not taken into account. There is another shortcoming of the 
simple expressions (JT|)-(P|: they do not describe the shift of the phonon frequency due to 
the interaction with disorder. Therefore, we have to use a more complicated technique of 
the Green functions jl^ . 

The heart of the technique is the phonon self-energy 

w \ _ v-^ Nk-q)| 2 , s 

The real functions and are now themselves defined by the real and imaginary 

parts of the self-energy. We obtain them by solving the system of the equations 

Q 2 {u) -u 2 - iu[T(u) - T nat ] - c d v{q = 0) = (6) 

The last term (named the homogenous shift) in the left-hand side of Eq. © gives the phonon 
shift due to the averaged effect of disorder, whereas the self-energy, Eq. (jSJ), describes the 
fluctuation effect of the defects. In addition to the phonon shift (inhomogeneous, depending 
on the frequency), the self-energy produces a phonon width. Note, that if the phonon does 
not interact with disorder (v = 0), Eqs. (jEJ) naturally give Q(u) = u and T(u) = r nat . 
Next, in the Born approximation, Eqs. ©-© give the simple expression, Eq. for the 
phonon damping. 

In the Raman light scattering, the wave- vector transfer q ~ lu^/c is determined by the 
wave vector of incident light, being the incident light frequency. In the integral ©, the 



values of k ~ l/l = ^Tluq/s are important, where / plays the role of the phonon mean free 
path. The dispersion parameter s~5x 10 5 cm/s is typically of the order of the sound 
velocity and the width T is of the order of (1CT 1 -j- lCT 2 )^. Thus, the condition q k holds 
and we can omit q in Eq. 



To simplify the calculations, we assume that the potential function v(k) takes constant 
value vq in the region k < 1/ro and is equal to zero for k > I/tq. Then the parameter ro 
gives the size of the region in the real space, where phonons are scattered by defects. If we 
consider the phonon scattering by strain fluctuations, r has the meaning of the correlation 
length of these fluctuations. 

The line asymmetry is very sensitive to the relation between the parameters tq and I. 
The case where ro < I is referred to as the short-range disorder. It is easy to verify that the 
line asymmetry is much larger for the short-range disorder than in the opposite case. We 
can rewrite the condition of the short-range disorder as 



IT < 1, 

using the estimate of the Debye frequency uoo = its /a, where a is the lattice constant. 
For the point defects, the calculation of £(a>), Eq. (0), gives 

" . . ,/l (b + a^ + a 2 

b- (ai -ra 2 )(-ln— — — ^ 

v 4 (o — ai) 2 + a 2 



(7) 



E(u) = A 



(8) 
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in the case of maximum of the phonon branch and 



T,(u) = A 



, , /l (b + a 2 ) 2 + aj 
-b + f o 2 + iai)\ - In — k 

v ' V4 (h — a„\ 2 4- n 2 



% b + a 2 

-- arctan 

2 ai 



2 i „2 
1 

(b - a 2 ) 2 + a\ 

% b - a 2 

- arctan 

2 ai 



(9) 



in the case of minimum. Here 



ai = [tf(u)-u 2 + p(u)]V 2 /V2, 
a 2 = l-n 2 (u) +u 2 + p(u)] 1/2 /V2, 

pH = [(fi 2 M-^ 2 ) 2 + ^ 2 r 2 M] 1/2 , 

b = s/r a , A = c d v 2 /27r 2 s 3 . 

Notice, that the dimensions are the following: s [cm u], Cd [1/cm 3 ], w [cm 3 w 2 ], A [u]. 

Solving the system of Eqs. we find the functions fl(u) and T(u). Then we can 

calculate (see, i.e., the Raman intensity 

ujT(uj) 



[1 - exp(-huj/k B T)}[(n 2 (uj) - uj 2 ) 2 + uj 2 Y 2 (u) 
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Equation (fTUj) can be applied to the Stokes lines {u > 0), as well as to the anti-Stokes ones 
(u < 0). 

If = uq and T(u) = r nat are constant, Eq. (|l(Jj) gives the Lorentzian line-shape, 

because T nat <C ujq. But, as we have already seen, the non-Lorentzian line-shape is obtained 
if the width T(u) depends on to. 
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FIG. 1: Raman spectra (points - the experiment data, solid lines - the theory) from the sample 
in as-deposited conditions (bottom) and after annealing for 989 hrs at 1200 °C (top). 

Up to this point we have considered a single phonon line. In the general case there are 
several lines with their positions and widths I\(u;). The corresponding terms have 

to be summarized in Eq. (fTTHl . The intervalley terms with Vij could be written also in the 
self-energy These terms are essential when the disorder removes the degeneration of 
branches. 



III. COMPARISON OF THE THEORY WITH THE EXPERIMENTAL DATA 
FOR ZIRCONIA 

Two Raman spectra taken from paper are shown in Fig. 1. They were recorded 
from two samples of yttria-stabilized (8.6 mol% YO1.5) zirconia - the first in as-deposited 
condition (bottom) and the second (top) annealed at 1200 °C. Narrowing and decrease of 
asymmetry of the lines are evident. The fit to the theory, Eqs. (JEI)-©) is shown in solid lines; 
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FIG. 2: Raman line 470 cm -1 (points - the experiment data, solid lines - the theory) from the 
sample in as-deposited conditions (left panel, bottom) and after annealing for 989 hrs at 1200 °C 
(left panel, top); the peak of line from the as-deposited sample is marked by the vertical dash- 
dotted lines. The asymmetry of the line indicates that the phonon branch has a minimum. In 
the right panel, the total width T(iv) (top) and the inhomogeneous shift Q(uj) — ujq (bottom) are 
plotted versus the frequency transfer for two spectra. 
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FIG. 3: Same as Fig. 2, but for the line 250 cm . The phonon branch now has a maximum. 
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parameters of the fit are listed in Table. In fitting, we assume that the effect of disorder on 
the spectra from the sample after annealing is negligible (A = 0). There is no information 
about the width of Raman lines for pure zirconia. Therefore, the line widths for the sample 
after annealing are referred to as the natural widths r nat . 

The broadening T — r nat of the line for the sample in as-deposited conditions is induced 
by the phonon-disorder interaction determined by the parameter A. We assume, that the 
natural width r nat of each line before annealing is equal to r nat of the same line after 
annealing (see Table). 

The line shape of spectra for the sample with disorder is determined by two factors: (i) 
the value of the parameter r and (ii) the extremum type (minimum or maximum) of the 
phonon branch. In Figs. 2 and 3 the results of fitting are shown in detail. The 470.5 cm -1 
line is shown in Fig. 2, left panel. The high-frequency wing of the line from the as-deposited 
sample is more gentle. This corresponds to the minimum of the phonon branch. In the 
right panel, the width T(u) and the inhomogeneous shift (Q(u) — uj ) are plotted on the 
top and bottom, respectively, as functions of the frequency transfer u. For the sample after 
annealing, they are constant = T nat and = ujq. 

We see that the width T(uj) of the line from as-deposited sample increases with the 
frequency transfer around the line peak Qq, which is marked by the dash-dotted vertical 
line. This resembles the square-root dependence in Eq. (j2J) for the minimum of the phonon 
branch and we use Eq. Q in fitting. At the peak position, the width takes the value r ; 
this value is given in Table, it is slightly less than the FWHM (full- width at half-maximum) 
because of the frequency dependent T(u). 

The behavior of the 250 cm -1 line (see Fig. 3) from the as-deposited sample has the 
opposite character: the low-frequency wing drops more slowly. Then we use Eq. (jHj) corre- 
sponding to the maximum of the phonon branch. The asymmetry parameter As listed in 
Table is a difference between the right and left wings at the half-maximum. This parameter 
takes positive (negative) values for the phonon-branch minimum (maximum). The asymme- 
try parameter and the FWHM were determined for the line extracted from the total spectra 
with the help of the fit. 

Several points can be noted in Table. First, the radius of defects (correlation radius) 
is about the value of the lattice parameter a. This is the case of a short-range disorder 
resulting in the evident asymmetry of the line shape. Second, because of disorder, the line 
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TABLE I: Final values obtained in fitting for spectra from the s amp les after annealing at 1200 °C 
(Ann) and as deposited (Asd): the experimental (according to [l3|]) phonon frequency Q exp , the 
peak position 0,q, the peak position uq including only the homogeneous shift, the total line width 
To (subscript is referred to as the peak position), the natural line width r nat , the line asymmetry 
As, the interaction parameter A (all in cm -1 ), the defect radius ro/a (dimensionless), the relative 
intensity Iq, and the type of extremum. 



line 


^exp 


O 


UJ 


r 


pnat 


As 


A 


ro/a 


h 


extr 


1 Ann 


149 


148.6 


148.6 


8.7 


8.7 








- 


0.215 


min 


1 Asd 




151.2 


152.0 


14.6 


8.7 


3.1 


225 


1.60 


0.19 


min 


2 Ann 


269 


264.0 


264.0 


25.5 


25.5 










3.32 


max 


2 Asd 




250.9 


245.0 


48.0 


25.5 


-23.7 


225 


0.73 


2.55 


max 


3 Ann 


319 


321.0 


321.0 


26.0 


26.0 










1.15 


max 


3 Asd 




330.4 


325.0 


45.0 


26.0 


-18.4 


200 


0.83 


1.10 


max 


4 Ann 


461 


463.3 


463.3 


20.0 


20.0 










1.55 


min 


4 Asd 




470.5 


472.5 


26.5 


20.0 


5.2 


100 


1.27 


1.45 


min 


5 Ann 


602 


604.5 


604.5 


36.0 


36.0 










3.30 


min 


5 Asd 




602.1 


605.0 


43.3 


36.0 


5.8 


75 


0.95 


2.30 


min 


6 Ann 


648 


643.0 


643.0 


22.0 


22.0 










6.70 


max 


6 Asd 




636.4 


633.5 


27.2 


22.0 


-3.9 


50 


0.95 


7.85 


max 



position is slightly shifted (by less than 5 %), but the line width increases considerably. This 
can be understood with the help of Figs. 2 and 3 (right panels). The frequency dependence 
of the shift function Q(u) — u has a zero about the line peak, whereas the width function 
r(u;) takes almost its maximum value. 

IV. DISCUSSION 

Lughi and Clarke explained their results in the following way. The Raman line becomes 
more symmetric after annealing. According to the condition ((Zj), this means that the cor- 
relation radius of disorder in as-deposited samples must be quite small, of a few lattice 
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parameters. Therefore, the disorder can not be induced by the large-scale strain fluctua- 
tions, for instance, of the crystallite size (about 50 nm). 

The disorder with small correlation length can be realized by Y 3+ ions and their associated 
oxygen vacancies. But the concentration of yttria (8.6 mol%) before and after annealing is 
the same in the sample. Lughi and Clarke proposed a redistribution mechanism of the Y 3+ 
ions and oxygen vacancies. They reasoned that the regions of the c phase become richer 
(and the t phase is more pure) in the Y 3+ ions and oxygen vacancies in the t — c phase 
transformation during annealing. Because of this redistribution of defects the Raman active 
t phase exhibits the symmetric and narrowing Raman lines. 

Considering the ions and vacancies as the point defects that scatter the long-wave optical 
phonons, we use Eqs. (0)-© in our comparison. 

Finally, as can be seen from Table, the phonon-disorder interaction is larger for the low 
frequency lines. This interaction decreases progressively with the frequency of the lines. We 
conclude that the distribution of the heavier defect ions, i.e., Y 3+ , affects the phonon modes 
stronger than the disorder of oxygen vacancies. 
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